On the Capture of Dark Matter by Neutron Stars 
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We calculate the number of dark matter particles that a neutron star accumulates over its lifetime 
as it rotates around the center of the Galaxy, when the dark matter particle is a self-interacting 
boson but does not self- annihilate. We take into account dark matter interactions with baryonic 
matter. We show that dark matter self-interactions play an important role in the rapid accumulation 
of dark matter in the core of the neutron star and its collapse. For the case when the dark matter 
density is 1 GeV/cm^ and the dark matter self-interaction cross section is 10~^''cm^, the observation 
of a typical neutron star implies strong limits on the dark matter-nucleon cross section for a broad 
range of dark matter masses. We determine the excluded region of the parameter space for dark 
matter mass and dark matter interaction cross sections based on the observation of old neutron 
stars. We discuss the dependence on dark matter density. We show that for a dark matter density 
of 10'^ GeV/cm^, the constraint on the dark mattter self- interaction cross section is several orders 
of magnitude stronger than the current Bullet Cluster limit. 
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I. INTRODUCTION 

There is overwhelming evidence of the dark matter ex- 
istence in the Universe, from the observation of miss- 
ing mass in galaxy clusters [l[ to the precise measure- 
ments of the cosmological baryonic fraction performed by 
WMAP0I and BOSSfl. The possibility that the stan- 
dard gravitation law needs to be modified to explain the 
observations with the ordinary visible baryonic matter 
has recently been ruled out by the Bullet Cluster data 
[1]. The particle physics interpretation of dark matter 
requires dark matter particles to be weakly interacting 
and in thermal equilibrium until the Universe expansion 
becomes such that particles cannot find each other and 
their interactions freeze-out. Large-scale structure for- 
mations indicate that dark matter particles need to be 
non-relativistic at the time of the freeze-out, i.e. dark 
matter needs to be "cold." Measurement of the matter 
density and its baryonic component implies that dark 
matter density contribution is about 25% [l-Ql. 

Since the dark matter density is inversely proportional 
to the dark matter annihilation cross section at freeze- 
out, the observed density of dark matter in the Universe 
today constrains the annihilation cross section in the 
early Universe, specifically at the time of freeze-out. On 
dimensional grounds, a dark matter particle with mass 
in the range of 100 GeV to several TeV with weak scale 
couplings can have annihilation cross sections of the or- 
der of {av)ann = 3 X 10~'^^ cm'^ at the freeze-out, 
providing a natural explanation for the observed density 
of dark matter today There have been recent dis- 
cussions of the possibility of asymmetric dark matter, 
where the dark matter particles are not self-conjugate, 
see for example, Refs. j8l-[l2| and references therein. An 
initial particle-antiparticle asymmetry ultimately leaves 
non-annihilating DM particles remaining in the current 
epoch. This is the case we consider in this paper. 



As an astronomical object in the Galaxy rotates 
around the center in its orbit, it will sweep through the 
Galactic dark matter halo and eventually capture some 
of the particles on its way. In time, dark matter parti- 
cles that are captured may have effects on the observa- 
tional properties of the astronomical object, which may 
then be used to constrain the nature of the dark matter 
pT| - [23j . In that respect, neutron stars provide a natu- 
ral laboratory to constrain the properties of dark matter 
[ill . [l3-im • Even though the surface area of a typical 
neutron star is much smaller compared to a more tradi- 
tional astronomical object, like the Sun, two properties 
make neutron stars very efficient in capturing Galactic 
dark matter particles. First, the immense baryonic den- 
sity inside a neutron star provides a natural area where 
there it is very likely that dark matter particles will in- 
teract and lose energy. Second, because of the strong 
gravitational force, it is also almost impossible for a dark 
matter particle to escape from a neutron star once it loses 
some of its energy through interactions. 

It is then only a matter of time for a neutron star 
to capture enough number of dark matter particles that 
will affect its observational properties. If the dark matter 
particles are annihilating, one such effect can be seen in 
the cooling of an old neutron star. The energy outcome 
of the annihilation process will result in an increase of 
the temperature that will remain constant and discern- 
able from other cooling processes in time (see, e.g., Ref. 
p^). Calculations of the annihilation effects on the cool- 
ing of a neutron star show that the resulting effective 
temperature of a neutron star would be approximately 
3000—10000 K [l^, depending on the local dark matter 
density, and the mass and the radius of the neutron star. 
However, the emission of a blackbody at these tempera- 
ture peaks at the UV to optical wavelengths, where the 
Galactic extinction hampers our observational capabili- 
ties to obtain precise measurements of the surface tem- 
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peraturcs of neutron stars unless they are very close. 

However, if the dark matter is not annihilating, under 
certain conditions, the capture process may have observ- 
able effects. For some values of the local dark matter 
density, dark matter mass and its interactions with nu- 
cleons and amongst themselves, the number of particles 
may reach to a value where the dark matter becomes self- 
gravitating. Once a dark matter inside a neutron star be- 
comes self-gravitating, and occupies small region in the 
core of the neutron star, it may collapse into a black hole, 
which would destroy the whole neutron star. In such a 
case even the very existence of neutron stars at certain 
ages can be used to constrain the properties of dark mat- 
ter. This has been studied in case when dark matter has 
interactions with baryonic matter only [ll|, [13, [l^ . We 
focus on the effect of dark matter self interactions. 

We start by reviewing the conditions necessary for a 
neutron star to collapse as it captures dark matter par- 
ticles. We discuss conditions necessary for dark mat- 
ter particles to be relativistic as they become captured 
by the neutron star, for the case when the dark mat- 
ter particle is a fermion or a boson. In section II we 
discuss thermalization of the captured dark matter par- 
ticles. For bosonic dark matter, we review conditions 
for self-gravitation and discuss the case when bosonic 
dark matter forms a Bose-Einstein condensate. In sec- 
tion III, we outline our calculation for the time evolu- 
tion of the dark matter particles as they get captured 
by a neutron star. We review experimental constraints 
on dark matter-nucleon and dark matter self-interaction 
cross sections. We discuss theoretical uncertainties in 
dark matter density profiles. In section IV we summa- 
rize our inputs and parameters used in the calculation. 
We give a simple example and rough estimate for the 
range of dark matter cross sections that would result in 
providing dominant effects on the increase of the number 
of captured dark matter particles as evolved in time. We 
discuss our results in section V, followed by conclusions 
in section VI. 



II. CONDITIONS FOR THE COLLAPSE OF 
NEUTRON STAR 



A. Chandrasekhar limit 

Our focus is on dark matter collapse to a black hole 
within a neutron star. As discussed below, thermaliza- 
tion of the dark matter particles is an important feature 
for fermionic dark matter. Subsequent accumulation to 
the limit of self-gravity is important if bosons have an 
impact on neutron star collapse. We begin by reviewing 
the well-known results for neutron star collapse. 

For neutron stars, collapse occurs only when the nucle- 
ons are relativistic. The relativistic energy of the neutron 



star is approximately 
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where G is Newton's constant and m„ is the neutron 
mass, and Nn is the number of neutrons. When the grav- 
itational energy of the neutrons equals the energy due to 
the relativistic Fermi momentum (E = 0), the number 
of neutrons is determined, independent of the radius r 
of the neutron star. Additional neutrons cause collapse 
into a black hole. The limit on Nn when E — is the 
Chandrasekhar limit for neutron stars, 
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More generally, for fermions with mass to^^, the Chan- 
drasekhar limit is 
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where the interactions with neutrons are neglected. An 
essential feature in this derivation is that the fermions are 
relativistic. The average energy of a relativistic fermion 
must satisfy 
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This translates to a requirement that the number of 



fermions Nf should be larger than A^p', 



N} > 1.6x10' 
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in order that the fermions are relativistic. The num- 
ber of neutrons in a neutron star with a typical radius 
of i? = 10.6 km and neutron star mass M = IAAAIq 
is such that the neutrons arc non-relativistic. For dark 
matter fermions to cause the neutron star to collapse, the 
requirement for relativistic energies may not be satisfied 
unless r <C 10.6 km or the dark matter mass is small 
compared to m„. The typical radius containing most of 
the dark matter reduces when dark matter is thermalized 
with neutrons in the neutron star, in which case r <^ R 
is the thermalization radius rather than the neutron star 
radius. We return to this below. 

For bosonic dark matter, there is no Fermi pressure. 
The corresponding Chandrasekhar limit for bosons, again 
neglecting the gravitational energy associated with neu- 
trons, is [ll| 



Ch 
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1.5 x 10 



34 



100 GeV y 



(6) 



since the kinetic energy per boson is of order E ^ hc/r, 
the zero point energy due to the uncertainty principle. 
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This approximate limit is confirmed by a number of ap- 
proaches to calculating the Chandrasekhar limit j24| . 
The introduction of self-interactions modifies this limit 
for bosons, potentially raising the limit to close to the 
Chandrasekhar limit for fermions [25j, JGj , although this 
is model dependent. For this paper, we rely on Eq. © 
for the Chandrasekhar limit for bosons. 

If the dark matter particles are bosons, they are not 
relativistic unless they get trapped inside the small re- 
gion of the neutron star. The requirement for relativistic 
bosons is that E ^ hc/r > mc^, so the relativistic condi- 
tion is 



100 GeV 10.6 km 



< 2 X 10-^ 



(7) 



After the thcrmalization time, bosonic dark matter is 
primarily within a small radius, but that radius is not 
small enough for Eq. ([7]) to be satisfied. A larger dark 
matter density is required for the bosons to form a black 
hole in the core of the neutron star. An additional stage 
of accumulation for bosons occurs when bosons form a 
Bose-Einstein condensate [lll[T7l . [20| . Finally, when p,^ > 
Ph, the dark matter is self-gravitating, namely, the baryon 
density pb can be neglected. During this final stage of 
collapse, the bosonic dark eventually becomes relativistic 
and may satisfy the Chandrasekhar limit. 

We now turn to thermalization and two scales, the 
thcrmalization radius in which the dark matter particles 
are found when they have thermalized, and the thermal- 
ization time, which should be shorter than the age of the 
neutron star to be relevant for dark matter collapse of 
neutron stars. We also review the scale of Bose-Einstein 
condensation for bosonic dark matter. 



B. Thermalization/condensation of the captured 
dark matter particles 

Thermalization occurs as dark matter particles inter- 
act with the nucleons inside the neutron star. Through 
interactions, they lose energy and eventually thermal- 
ize. The thermalization time depends on the dark matter 
mass, the temperature of the neutron star and the proba- 
bility of dark matter interaction with nucleons expressed 
in terms of the cross section 
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For light dark matter, < m„ jlll |. 
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tth = 7.7 X 10"^ yrs 



I GeV\ /2.1 X lO^'^'^cm^ 
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Here T is the temperature inside the neutron star. Once 
the dark matter particles reach thermal equilibrium they 
are contained in the core of the neutron star within the 
thermal radius given by [l7j 



9T 



Tth 



24 cm 



V 105 K 



100 GeV y/^ 



(10) 

where pb is the baryon density inside the neutron star, 
which we take tobel.4xl0^5 g/cm^, and tT^„ is the dark 
matter-nucleon interaction cross section. 

With thermalization, for fermionic dark matter, one 
can use rth in eq. ([5|) to find that 



Nf>Np'{r = rth) = 1.9 xlO^^] 



T 



3/2 



100 GeV lO^K 



(11) 

in order for fermionic dark matter to be relativistic. For 
T = lO^K, the Chandrasekhar limit for fermions, which 
requires relativistic particles, is applicable only for dark 
matter masses below 100 GeV when t > tth, since for 
< 100 GeV, Nf" > N^f^K When > 100 GeV, 



fermionic dark matter within the thermal radius with 
Nf ~ N^^ is non-relativistic. This appears not to have 
been taken into account in Ref. (22| . 

Fig. [1] shows the relevant Nf for the Chandrasekhar 
limit, for N'^^^^{r = rth) £md the number of fermions re- 
quired for relativistic fermions when r = rjy, the radius 
of the neutron star. Before the thermalization time, the 
larger of N^^ and N'j'^\r = r^) determines whether or 
not the dark matter fermions collapse to a black hole. For 
t > tth, it is the maximum of Nj"^ and N'f'^\r = rth)- 

If the dark matter particles are bosons, there are ad- 
ditional stages of accumulation of dark matter parti- 
cles. Within the thermalization radius, as the number of 
bosons continues to increase in time, the bosons can fur- 
ther form a Bose-Einstein condensate (BEG) [l7|. The 
particle number density required to form a BEG in a 
sphere of radius rth in the neutron star is p7l [20| 
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(r = rth) = 2.7 x 10^ 
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105 K 
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Once the BEC forms, the radius in which the dark mat- 
ter particle can reside, which is the radius of the wave 
function of the ground state in the gravitational poten- 
tial of the neutron star becomes even smaller than rth 
[Ulilli^, namely. 
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100 GeV 
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(13) 



Clearly, rsEC < fth, which has significant implications 
for reaching the Chadrasekhar limit in less time. 
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FIG. 1: As a function of fermionic dark matter m^, the logio 
of the number of fermions for the Chandrasekhar Hmit (which 
requires the dark matter to be relativistic) shown with the 
solid line and the minimum number of fermions required for 
relativistic energies when r = rjv (dot-dashed) and r = rth 
(dashed). 
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FIG. 2: As a function of bosonic dark matter m^, the logio of 
the number of fermions for the Chandrasekhar limit (which 
requires the dark matter to be relativistic) shown with the 
solid line, the minimum number of fermions required for the 
bosons to be self-gravitating when r = rth (dot-dashed) and 
r = tbec (dotted), and the number of bosons within the ra- 
dius r = rth for a Bose-Einstein condensate to form (dashed). 



C. Conditions for bosonic dark matter 
self-gravitation 

If the dark matter particles are bosons, they can only 
become relativistic and collapse if they become self- 
gravitating. This happens when the density of the dark 
matter particles exceeds the baryon density within the 
same volume. Since the nucleon density does not change 
in time in the neutron star, as soon as the number of dark 
matter particles reaches a critical number, the bosons be- 
come self-gravitating. This critical number is given by 



self 



For r = rth, 
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1 GeVy/^ 
) 
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105 K 



(14) 



3/2 



(15) 

Clearly, as the radius of the dark matter core decreases, 
the condition for the dark matter particles to become 
self-gravitating will be easier to reach, which would even- 
tually cause the dark matter core to collapse and form 
a black hole. When the Bose-Einstcin condensate has 
formed. 



Nselfir = TBEC) - 1-1 X 10 



28 



1 GeV 



5/2 



(16) 



Comparing the equations for N^^ (eq. ^) and 
for self-gravitation in eq. (fT51) one can see that the 
Chandrasekhar limit is already exceeded for thermalized 



bosonic dark matter when self-gravity is established, as 
shown in Fig. [2l Self-gravity results in further infall and 
ultimately relativistic bosons, the final requirement for 
Eq. dS]) to apply. 

To summarize, the necessary condition for collapse of 
neutron star to a black hole is that the dark matter is 
thermalized within the age of the neutron star, which 
we take to be 10^ years. Additional condition that has 
to be satisfied in order to ensure the collapse is that 
the number of captured dark matter particles exceeds 
following limits: N^^ for < 7.4 GeV, N^^^ for 
7.4GeV < < 1.26 x 10^ GeV, and Nseif{r = nh) for 
> 1.26 X lO'* GeV. Clearly, if the dark matter par- 
ticles are bosons, the collapse into a black hole will be 
much easier than for fermions. Once a body of mass be- 
comes self gravitating and there is no force to act against 
gravity, its outcome is determined. 



III. CAPTURE OF DARK MATTER 
PARTICLES BY NEUTRON STARS 

The time evolution of the dark matter particles cap- 
tured by a neutron star is given by [l^ | 



dt 



= Cc + CsN^-CaNl 



(17) 



where Cc is the capture rate due to dark matter-nucleon 
interactions, CgN^ is the capture rate due to the dark 
matter self-interactions and CaN^ governs the number 
of particles lost due to the their annihilation. We con- 
sider the case of asymmetric dark matter, when dark 
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matter particles do not annihilate. In this case, for time- 
independent Cc and Cs, the solution for the Eq. (17) 



1) 



(18) 



A. Dark matter-nucleon interactions 

The number of dark matter particles that can be cap- 
tured by a neutron star is given by [T9| 



^ 2 Px f 3 
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(19) 



where the general relativistic affects assuming a 
Schwarzchild geometry have been incorporated. Here 
is the dark matter density at the location of the neutron 
star, is the mass of the dark matter particles, M and 
R are the mass and the radius of the neutron star, re- 
spectively, and V is the average velocity of dark matter 
particles in the Galactic halo [l^. For a dark matter 
particle to be trapped with a single collision, E < niy^eo, 

which defines eq- Since eo ^ we take e"^^"/'" ~ 
p^ . In order for a dark matter particle to stay inside 
a neutron star it must lose enough energy by interacting 
with the particles inside the neutron star. The last fac- 
tor in eq. (fT9)) . /, determines the fraction of the dark 
matter particles that will be trapped inside the neu- 
tron star due to interactions with nucleons, which for 
< CTcrit = 7rm„i?^/M ~ 2 X 10~^^ cm^ can be ap- 



proximated as 

/ = (l - cxp 



-dl 



-dl 



(20) 



where dl is the infinitesimal arc length of a trajectory 
in the neutron star [isl . [l9| . Assuming constant density 
for the neutron star, numerical evaluation of / is terms 
of the critical cross section, iJcrit, which is geometrical 
cross section of the neutron star divided by the number 
of nucleon targets, is given by [l^ 



/ 



'Xn 



dl 



0.45- 



^crit 



(21) 



This equation is applicable to the case where (t^„ < acnt , 
however / saturates to unity when (t^„ is larger than a„it 
p^ . Thus, the quantity / is 



/ ~ 0.45 



min((T;^„, CTcrit) 



(22) 



Currently, the most stringent experimental limits on 
the DM-nucleon cross section, a^n, come from the 
XENON [13 and CDMS [H experiments which look for 
energy deposition via nuclear recoils from dark matter 



scattering. For dark matter mass between 20 GeV and 
100 GeV, the limit is a^n < 6x10^''^^ cm^, while for larger 
mass, i.e. = 10'^ GeV, the limit is about an order 
of magnitude larger and similar for lighter dark matter 
particle, when ^ 15 GeV. Recently, CRESST-II has 
claimed A.7a signal corresponding to niy, ~ 10 ^ 30 GeV 



and for 10 



-40 cm2 < a. 



< 3 X 10-*^cm^ [291. This sig- 



nal docs not seem to be consistent with DAMA results 
claiming annual modulation effect consistent with dark 
matter particle of mass around 10 GeV with the cross 
section a^n ^ 10^*°cm^ [s^j- We will incorporate these 
limits by taking conservative approach and considering 
only cases when a-^n is below lO""** cm^. 



Name 



Value 



Solar Mass 
Velocity dispersion of DM 

Distance of the NS 
from the Galactic Center 
Neutron Star Mass 
Neutron Star Radius 
Central Density of Neutron Star 
Temperature inside the Neutron Star 



Mq = 1.98892 X lO'^'^g 
V = 220 km/s 
d = 3.7 kpc 



M = 1.44 M© 
R = 10.6 km 
Pb = 1.4 X lO^^g/cm' 
10^ K 



TABLE I: Constants and parameters used in the results pre- 
sented here. 

The parameters that we use are given in Table HI With 
these parameters, the capture rate, Cc, for a^n < f^crit = 
2 X 10""*^ cm^, is given by. 



Cc = 9.19 X 10-^-^ 
m^/GeV 



(23) 



B. Dark Matter Self Interactions 

Once a certain amount of dark matter particles accu- 
mulates in the neutron star, their very existence inside 
the neutron star will affect the capture of new dark mat- 
ter particles due to their self- interaction [l3|. The self- 
capture rate, Cs, is given by [l^ 



Cs 



\zZ-'^XxVesc{R) 



,{R),i ^erf(r;) 



1 



2 m 



V 1 



2GM 



(24) 

where, cr^^ is the dark matter elastic scattering cross- 
section and Vesc{R) is the escape velocity from the sur- 
face of the neutron star. For a 1.4M0 mass neutron star 
with R = 10.6 km, the escape velocity from the surface is 
Vesc{,R) = 1.9 X 10^° cm/s. We modified eq. ^ follow- 
ing Ref. [l^ in order to take into account the general rel- 
ativistic effects assuming as Schwarzchild geometry. The 
quantity (j)^ which signifies how c omp act the star, is a 
dimensionless potential defined as [Tj], 



escV } 



(25) 
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We take {4>^) = 1, which is valid in the approximation 
that the mass density of the neutron star is uniform. 
Since the density is larger in the core, which implies 
that {(j)^) > 1, our assumption is conservative. Finally, 
77^ = 3/2{vn /v)'^ depends on the velocity of the neutron 
star vn in the Galaxy. We approximate erf (77)777 ~ 1. 
With our choice of the parameters given in Table I, Cg is 
given by, 



TTix/GeV 

Currently the most stringent limits on self-interaction 
cross section, cr^^, come from the Bullet Cluster obser- 



vations I31|, i.e. 



< 2 



(27) 



Dark matter density profiles are usually 

parametrized as a function of the distance r from 
the Galactic Center. 



(28) 



where ps = 0.26 GeV/cm'^ and the values for the param- 
eters a, (3, 7 and rs for several dark matter profiles are 
given in the Table II. 



Model 


a 


P 


7 r 


■s (kpc) 


Navarro-Frenk- White (NFW) [32] 


1 


3 


1 


20 


Moore [33] 


1.5 


3 


1.5 


28 


Kravtsov [34] 


2 


3 


0.4 


10 


Isothermal with core radius \35} 


2 


2 





3.5 



TABLE II: The values of the parameters in Eq. for 
different DM density profiles. 

In our calculation we use Navarro-Frenk- White (NFW) 
dark matter halo profile [s^j- At a distance of d = 3.7 
kpc from the Galactic Center, the dark matter density 
according to the NFW halo profile is p-^ = 1 GeV/cm'^. 
The uncertainty due the choice of the dark matter den- 
sity profiles depends on the location of the neutron star 
relative to the Galactic Center. For example, for distance 
of r = 5 kpc, dark matter density obtained from Moore 
[33! would give about 10% higher value than NFW, while 
Kravtsov [3^ would be about 25% lower than NFW. This 
uncertainty is larger for neutron stars closer to the Galac- 
tic Center. 



IV. DARK MATTER CAPTURED BY THE 
NEUTRON STAR 

We start by considering the condition for which the 
self-interactions will play a role in the neutron star col- 
lapse. To estimate the time at which dark matter self- 
interactions become important, we impose the condition 



that Cst ^ 1. At this time, given by Eq. (18) will 
change its time dependence from being linear, — Cc t, 
to becoming an exponentially increasing function of time. 
Using the expression for Cg given by Eq. (28), we get. 



t(yr) ~ 0.94 x 10 



-2i "^x/GcV 
/cm^ 



(29) 



'XX 



If this time is reached for t ^ ttt, which is satisfied 



(26) when {my.GeV) / {a^^/ cnr) < 10^° for tth = 10 years 



then the exponential increase of will continue until 
the dark matter self capture reaches geometrical limit, 
i.e. a^y,N^ = ttR^. The time at which the geometrical 
limit is reached can be obtained by solving. 



Using the expression for Cc and Cg we get 
t f r) - ^ \n'^'^'^ ^ ^'^'^^ 



(30) 



(31) 



After has reached this geometrical limit, its growth 
with time changes to a linear behavior given by. 



= {Cc + Cs N^)t 



where 



CsN^ (yr-i) = 



33 



3.74 X 10 

m^/GeV 



(32) 



(33) 



and Cc is given by Eq. (23). This linear increase con- 
tinues until the thermalization time, at which point the 
linear increase in time is slowed down due to the several 
orders of magnitude decrease in CgN^, which is now ob- 
tained with the condition that a^^N'^ = '^''"'th- 

Fig. [3] illustrates the impact that dark matter self in- 
teractions have on the time evolution of the number of 
captured dark matter particles. As an example, we have 



taken 771-1 



10 GeV, fj^ 



10 2^ cm^ and a^n = 10 



■55 



cm . Using Eq. (8) we find that thermalization time is 
10^ years. Initially the number of captured dark matter 
particles increases linearly with time due to dark matter 
interactions with nucleons. At i = 0.94 x 10^ years, dark 
matter self-interactions become important. This can be 
seen by solving Eq. (29) for the parameters chosen in 
the figure. This exponential increase of ceases when 
the geometrical limit for cr^^ is reached. Using Eq. (31) 
above, we find that this happens at to = 2.3 x 10^ years. 
From this point on, the number of captured dark mat- 
ter particles has a linear increase with time, according to 
Eq. (32). When t = 10^ years, dark matter thermalizes 
within r = rth region. Since the number of particles at 
this point in time is larger than N^^'~^{r = rth), there is a 
further collapse of bosons to a Bose- Einstein condensate. 
The accumulation of dark matter particles is very slow at 
this time since there is a very small cross sectional area 



7 




10' 10^ lO' 10' 

Time (years) 

FIG. 3: The number of dark matter particles accumulated in 
a neutron star as a function of time, for ~ 10 GeV and 



= 10" 



cm , neglecting dark matter self interactions 



(lower green line) and with a^^ 
line). 



10" 



cm (upper blue 



of the neutron star intersected by the ambient dark mat- 
ter density. For t > Uh, -/V^ rises very slowly with time, 
because the increase is limited by the geometrical cross 
section for r = tbec- Note, however, that for t > ttjn 
> N^^ and > Ngei / , so bosonic dark matter with 
these interactions and mass will collapse to a black hole. 



V. RESULTS AND DISCUSSION 

For our choice of the parameters, the time to reach 
geometrical limit, tQ, defined by Eq. (30) and the ther- 
malization time, tth defined by Eqs. (8) and (9) are ap- 
proximately given by. 



= 0.9434x 10- !%Z^, 1,(1:^141^ (34) 



and 



= 1.136 X 10-5° rn^VGcYj ^ ^ ^^^^ 



o-^ri/cm^ 



(35) 



it/, = 1.6x10-5" ^— ] — - formx<lGeV. 

m^/GeV cr^^n/cm"^ 

(36) 

With these expressions and Eq. (18), the number of cap- 
tured dark matter particles can be obtained analytically. 

The number of captured dark matter particles at time 



t, for the case when tth > tQ, is given by 



9.19 X 10 



77 



CF^n/CTC? 3.74 X 10^3 



m^/GeV 



X (t-tG)/yr + 8.67 X 10 



56 



m^/GeV 
cr^n /cm^ 



X / 1-06x10^^ Wen.- Wyr _,p,) 
V mx/GeV ) 

Using the expression for tg, the Eq. (37) can be simpli- 
fied, 

/ N r„ „ ^770-yn/cm2 3.74 xlO^^,, 
N^{t) = [9.19 X lO^^^^^^f— - + ^^](t-tG)/yv 



TO^/GeV mJGeV 



12 



3.53 X 10 

crxx/cni2 



8.67 X 10 



56 "X" 



(38) 



When tth < to, the number of captured dark matter 
particles is given by, 



N^{t) = 8.67 X 10^ 



[expl 



1.06 X 10^ 



m^/GcV 



lp9) 



The quantity is evolved only up to the thcrmalization 
time because the increase in number of captured particles 
beyond tth to the final time of 10^ years is negligible due 
to the very small thcrmalizcd region. To determine what 
values of the dark matter cross section and its mass are 
excluded by the observation of old neutron stars, we im- 
pose the conditions for the collapse described in Section 
II. 

In Figures 01 [5] and [HI wc show the time evolution of the 
number of dark matter particles captured by a neutron 
star for = 1 GeV, — 10 GeV and = 100 
GeV, and for several values of dark matter interaction 
cross sections. We have taken the Galactic dark matter 
density to be 1 GeV cm"^ and other parameters as given 
in Table [H The self-interaction cross section is taken to 
be cr^^ = 10~2^ cm^ and wc consider two different values 
for dark matter nucleon interaction cross section, cr^„ = 
10-55 cm^ (left panels of Figs. 4-6) and cr^n = lO-''^ 
cm^ (right panels of Figs. 4-6). When the dark matter- 
nucleon cross section is small, (t^„ = lO'^s cm^, the dark 
matter self-capture is important, as shown by the upper 
solid curve labeled "C^" in the left panels of Figs. 4, 5 
and 6. Neglecting self-capture yields the lower solid line 
labeled "Cc." In the right panels of Figs. 4-6, cr^„ = 
10""'* cm^ and the self-capture contribution is negligible 
compared to the capture by nucleons, shown by the solid 
line labeled "Cc." 

Fig. d] shows our results for = 1 GeV. In addition 
to the solid curves showing the number of dark matter 
particles as a function of time, we also show the var- 
ious limits on the number of fermions and bosons dis- 
cussed in Sec. II. The uppermost (light blue) dashed line 
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represents N'j^^. For = 1 GeV, for our parameter 
choices, both panels show that < N'^^ , so fermionic 
dark matter is unconstrained for these parameters. The 
middle (red) dashed line represents N^^. In both pan- 
els, at t ^ 10^ years, > N^''', and TV^ > Nseif, 
so bosons become relativistic and the collapse occurs for 
these parameters for bosonic dark matter. In the right 
panel of Fig. SI where (t^„ = 10~^^ cm^, the collapse 
is entirely due to the capture of dark matter by neu- 
trons. The dashed (orange) line corresponds to Nseif- 
The value of Nseif decreases at the thermalization time 
so that N-^ > Nsei f , so bosonic dark matter will continue 
the infall to become relativistic. In addition, around this 
time, a Bose Einstein condensate will form, further com- 
pressing the bosonic dark matter. The quantity Nbec is 
shown by the dot-dashed line. Finally, the Hawking radi- 
ation limit, which corresponds to the number of particles 
needed to form a black hole which will not evaporate due 
to Hawking radiation, is shown with the dotted line la- 
beled by N'g^ [lO]. The right panel in Fig. H shows that 
> N^''' > Ng'}} for t = 10^ yr, so for these param- 
eters bosonic dark matter will collapse to a black hole 
in the core of the neutron star, followed by neutron star 
collapse. 

In Fig. 5, we show our results for the time evolution 
of the captured dark matter particles for ~ 10 GcV. 
We note that in this case, for cr^^ = 10"^'^ cm^, and 
a^n = 10~^^ cm^, at the time when Chandrasekhar limit 
is reached, around 10^ years, dark matter density had not 
reached the value necessary for it to start self-gravitate. 
At about 10^ years, thermalization happens and dark 
matter particles start to self-gravitate, i.e., Ny, > Nseif 
causing the neutron star to collapse. If we increase dark 
matter- nucleon cross section to 10""'^ cm^, dark mat- 
ter capture becomes driven by the dark matter-nucleon 
interactions, and the Chandrasekhar limit is reached in 
about 10^ years. 

For the case when the dark matter mass is 100 GeV and 
the dark matter nucleon cross section is 10~^^ cm^, we 
note from Fig. 6 (left panel) that dark matter thermal- 
izes and starts to self-gravitate at a later time, around 10^ 
years, at which point the neutron star collapses. Increas- 
ing a^n to 10~'**cm^, makes dark matter self- interactions 
negligible relative to its interactions with neutrons. In 
this case, at the time of the thermalization, around 10^ 
years, the number of captured dark matter particles is too 
small to cause the neutron star to collapse. At i '-^ 10^ 
years, the Chandrashekar limit is reached, however the 
dark matter is not relativistic. Dark matter continues to 
accumulate, and they form a Bose-Einstein condensate 
around 5 x 10* years, at which point the number of dark 
matter particles exceeds N^^'~^ , causing the neutron star 
to collapse. 

To explore the full bosonic dark matter parameter 
space that is probed by the observation of the old neu- 
tron stars, we start by considering dark matter self- 
interaction cross-section values, a^^ = 10"^"* cm^ (left 
panel of Fig. [7]) and 10~^^ cm^ (right panel of Fig. 



[7]). For each value, we have determined the dark mat- 
ter mass and the dark matter-nucleon cross section a^„ 
for which the neutron star would collapse. The upper 
excluded region (red shaded area) in Fig. [7] corresponds 
to the case when dark matter capture by neutrons is 
a dominant process. For < 7.4 GeV, the lower 
limit on cr^„ vs comes from the requirement that 

> Nrmn = N^'"''\ For 7.4 GeV < < 1.26 x 10^ 
GeV, the lower limit on cr^„ vs comes from the re- 
quirement that A^^ > Nmin = NBEc{t = Uh), and for 

> 1.26 X 10^ GeV, the lower limit on CTp^„ vs comes 
from the requirement that A'^ > Nmin = Nseif{t = tth)- 
In all three cases, = Cct and we have taken t = 10^ 
years. We have also taken into account that the thermal- 
ization time has a different dependence for m^j, < 1 
GcV, as given by Eqs. (35) and (36). 

While we have not explicitly imposed the Hawking ra- 
diation limit given in Ref. [20|, we find that Nmin > 
TV^'^t = 6 X 10^6 (GcV/m^) for to^ < 4 x 10*^ GeV. In 

particular, for Nmin = N^^^ir = nh) = 2.7 x 10^*^ for 
7.4 GcV< < 1.26 x 10^ GcV, the Hawking radiation 
limit is satisfied. If one instead applied a limit from N^^ 
or Nseif (r = fth), the Hawking radiation limit would 
not be satisfied. In this mass regime, the higher Nmm is 
required to ultimately get a BEG, with self-gravitating 
bosons that become relativisitic. 

The lower excluded (red) region corresponds to addi- 
tional parameter space excluded due to the dark matter 
self interactions. The upper limit on on a^n vs m-^ comes 
from the same requirements on A^^ as in the case of dark 
matter captured by neutrons only, as described above. 
Two distinctive regions correspond to the case when 
N^itc) > Nrain for larger masses and N^{tth) > Nmm 
for smaller masses. As we decrease the mass, the up- 
per limit on a^n becomes mass independent reaching the 
value of 4 X 10~^^cm^. This can be seen from Eq. (37) 
for N^(tth) which has a dominant second term with an 
TO~^ dependence, the same dependence on dark matter 
mass as the Chandrasekhar limit, giving the mass inde- 
pendent value for (t^„ stated above. The limiting value 
a^n for very small is independent of cr-^,^. However, 
the upper value for for which this limit is reached is 
determined by a^^. 

The Bullet Cluster limit, Eq. (27), is also shown in 
Fig. [71 We note that the thermalization is crucial for the 
dark matter particles to become self-gravitating, thus the 
region of parameter space for which the thermalization 
does not happen cannot be excluded by considering the 
collapse of the neutron star. The (green) shaded region 
labeled "No Thermalization" shows the region of param- 
eters for which tth > 10^ yr. 



VI. CONCLUSION 

We have calculated the time evolution of the number 
of captured dark matter particles for different values of 
dark matter self-interaction and dark matter nucleon in- 
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10° lO' 10^ lO' lO' lO' lO' lO' 10° lO' 10'° lO" lO'^ 10° lO' 10^ lO' lO' 10^ lO' lO' lO' lO' 10'° lO" lO' 

Time (years) Time (years) 



FIG. 4: Time evolution of the number of dark matter particles accumulated in a neutron star, for = 1 GeV, ct^x = 
and a^ri = 10^^^ cm^ (left panel) and (7^„ — 10"^* cm^ (right panel) . In both panels, the green solid line corresponds to the 
dark matter capture without self-interaction. Blue solid line in the left panel shows the evolution of the number of particles 
when the self-interaction is taken into account. Chandrasekhar limit for bosons is shown with the dashed red line and the 
self-gravity limit with the dashed orange line. Dashed-dotted line show the number of particles needed to form a BEG, while the 
dotted line represents the number of particles needed to form a black hole which will not evaporate due to Hawking radiation. 




FIG. 5; Same as Fig. 4, but for = 10 GeV. 



teraction cross-sections. Assuming a NFW dark matter 
distribution, we have considered an environment where 
the dark matter density is 1 GeV/cm'^, which corresponds 
to a distance of about 3.7 kpc from the center of the 
Galaxy. We note that these assumptions on the dark 
matter density and/or distance as well as the the age are 
reasonable given the Galactic age and location distribu- 
tion. As an example, 5% of all the pulsars in the ATNF 
Pulsar catalog have a characteristic age that is greater 
than 10^ years. Similarly about 6% of all the pulsars in 
the same catalog arc located within the central 3.7 kpc 
in the Galaxy [37| . 

We have determined the parameter space where the 
number of dark matter particles captured by the neutron 
star that become self-gravitating (and thus relativistic) 



and exceeds the Chandrasekhar limit for bosons within 
10^ years. For bosons with > 7.4 GeV, the relativis- 
tic requirement is more stringent than the Chandrasekhar 
limit. In Figure [5] we show the regions that can be ex- 
cluded for a 100 GeV dark matter particle. The dark 
matter self-interaction cross sections larger than 10^^^ 
cm^ are excluded by the observations of the Bullet cluster 
[3l| . The parameter space where dark matter particles 
do not thermalize within the age of the neutron star cov- 
ers the region below (t^„ = 10~^^ cm^. The observations 
of neutron stars cannot provide any information about 
this particular region. 

The upper excluded region in Fig. [51 i.e., a^n > 10"'^^ 
cm^, is independent of the dark matter self- interaction 
cross section because in this region dark matter capture 
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FIG. 6: Same as Fig. 4, but for = 100 GeV. 
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FIG. 7: The region of the parameter space that is excluded by the observation of old neutron stars (red shaded areas) and by 
the Bullet Cluster observation (turquoise shaded area). The left panel corresponds to the case when g^x = 10~^^ cm^ while 
the right panel is for — 10~^^ cm^. The region of the parameter space for which thermalization does not occur within the 
age of neutron star is shown with the dark green shaded area. 



is entirely due to its interactions with nucleons. Taking 
into account the dark matter self-interactions allows us 
to find another region which can be excluded for this 
dark matter mass. This is the parameter space where 
the nucleon interaction cross section is between 10~^^ 
cm^ to 10^^"' cm^ and dark matter self-interaction cross 
section is greater than 5xl0~^^ cm^. As noted in Sec. 
V, this region can be excluded because for these values, 
thermalization takes a long enough time during which 
self-interactions plays a significant role in capturing dark 
matter particles. 

The results shown so far have been evaluated with the 
assumption that the dark matter density near the neu- 
tron star is = 1 GeV/cm'^, corresponding to a dis- 
tance between the neutron star and the galactic center 
d = 3.7 kpc using the NFW dark matter halo profile. 
For larger dark matter density p^, the constraints on the 



cross sections are enhanced with an approximate scaling 
with dark matter density. 

We begin with the lower exclusion regions. First, the 
thermalization lines do not depend on the dark matter 
density, rather they depend on features of the neutron 
star itself, so the green shaded regions ("No Thermal- 
ization") in Fig. 7 do not depend on p^. As already 
discussed, the red excluded region, for low dark matter 
nucleon cross sections and > 0.1 GeV come from 

Tlim 



Nxitc) > N^"' and tc < Uh- The upper edge of the 



red exclusion region for low dark matter-nucleon cross 
sections is where to = tth- In the right panel, the hori- 
zontal limit of the lower red exclusion region comes from 
N^itth) > N^"^. For the whole lower exclusion region, 
to a good approximation, the limits depend on (T^„/p, 
so if — 10 GcV/cm^, the lower red exclusion region 
increases by a factor of 10. 
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Dark Matter Self Interaction (cm^) 

FIG. 8: The parameter space that is excluded by the obser- 
vation of the old neutron stars (red shaded area) and by the 
Bullet Cluster observation (turquoise shaded area). The re- 
gion of the parameter space for which thermalization does not 
occur within the age of neutron star is shown with the dark 
green shaded area. 

The upper exclusion region comes from the require- 
ment that N^{t = 10^ yr) > Nmin, where the three 
separate slopes depend on whether the Chandrasekhar 
limit, the Bose-Einstein condensate number or the self 
interaction limit determines Nmm- The upper exclusion 
region is independent of cr^^, with N^(t) = C'ct. Since 
Cc oc p^a^n, when increases, the excluded region de- 
creases, for example by a factor of 10 when = 10 
GeV/cm3. 

When (T^x = 10"^'' cm^, for > 7.4 GeV, the gap 
between the two exclusion regions for the case of = 
1 GeW /cir? is about a factor of 1000 for a^^n- When 
Px ~ 30 GeV/cm'^, this window for the cross section is 



closed. In the right panel of Fig. 7, where a^^ = 10~^ 
cm^, the gap between upper and lower exclusion regions 
for high masses is on the order of 3 x 10^, so a density 
of Px ^ 2000 GeV/cm'^ would exclude the range of a^n 
values for which there is thermalization for > 7.4 
GeV. Similar scaling arguments can be made for the lower 
mass regions, which is dominated by the second term in 
Eq. (38), the term proportional to Cc, and thus has a 
linear dependence on p^. 

For stars in the core of a globular cluster, dark matter 
density could be as high as lO'^GeV/cm'^ [l^. In this 
case we can put a stringent limit on dark matter self- 
interaction cross section, < 10~^^cm^ for > 10~^ 
GeV, which is several orders of magnitude below the Bul- 
let Cluster limit. 

Finally we note that we have used Eq. ^ to constrain 
bosonic dark matter in neutron stars. A more sophis- 
ticated Chandrasekhar limit for bosons would include 
modification in the case of strongly interacting bosons, 
but the value of self-interaction cross section for which 
this happens is model dependent (25], [s^ . It would be 
interesting to explore specific models for strongly inter- 
acting dark matter which may be ruled out by the obser- 
vation of typical neutron stars. 
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